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Abstract. A dynamic process of repeating collisions of a pair of trapped neutral 
particles with weak spin-dependent interaction is designed and studied. Related 
theoretical derivation and numerical calculation have been performed to study the 
inherent coordinate-spin and momentum-spin correlation. Due to the repeating 
collisions the effect of the weak interaction can be accumulated and enlarged, and 
therefore can be eventually detected. Numerical results suggest that the Cr-Cr 
interaction, which has not yet been completely clear, could be thereby determined. 
The design can be in general used to determine various interactions among neutral 
atoms and molecules, in particular for the determination of very weak forces. 
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1. Introduction 

It is well known tliat tlie study of scattering is an important way to understand 
the interactions among particles. Experimentally, the incident particles with certain 
momentum would hit a target to initiate the scattering. The accuracy of experimental 
data depends on a precise control of bombarding energy. For charged incident particles, 
the initial momenta are usually imposed by electromagnetic forces which can be precisely 
tuned. For neutral incident particles, the initial momenta are in general difficult to 
control precisely. Furthermore, the interactions among neutral atoms or molecules are 
in general weak (say, the van der Waals force). Therefore, the precise determination of 
these interactions is difficult. Consequently, these interactions are usually determined 
in an indirect way (say, comparing related experimental data of spectroscopy and/or 
thermodynamics of a many-body system with theoretical results based on a model with 
a given set of parameters). In this way the associated theoretical calculations are usually 
complicated, and uncertainty might exist. 

On the other hand, the trapping of neutral atoms has been realized via optical 
traps since 1998 [H [2] . This techniques is in progress and fewer atoms can be trapped 
recently [3]. It might open a new way for studying the interactions. In this paper, 
an idea is proposed and related theoretical calculation is performed to show how the 
scattering with precisely controllable initial status can be realized in a trap. It turns 
out that, as we shall see, the collisions would occur repeatedly. Thereby the effect of 
each individual collision can be accumulated and therefore enlarged. This might lead 
to a better understanding of very weak interactions among neutral particles. 

Since the pioneer experiment by Greismaier, et. al. [4], the Bose-Einstein 
condensations of atoms with a larger spin (say, ^^Cr) become a hot topic. These 
condensates are a new kind of matter aggregation having the magnetic dipole-dipole 
interaction Vdd more than twenty times stronger than that of the alkalis family. 



where the strength Cd = ^Qii\g% / {An) with /xq being the magnetic permeability of 
vacuum, ixb the Bohr magneton, and gp the Lande g factor, Fj the operator of the 
spin of the i-th atom, and r = r2 — ri. Consequently, the spatial and spin degrees 
of freedom are coupled so that the conversion of spin angular momentum into orbital 
angular momentum can be realized. Thereby new physical phenomena (say, rotonlike 
behavior) might appear O |6]. In addition to the long range Vdd-, the short range 
interaction is also spin-dependent and, for low-energy systems, can be in general written 
as = 5{ri — ^2) Xlsfi's^P"^' where S is the total spin, gs is the strength related to 
the s-wave scattering length of the S'-spin channel, and the projector of the S- 
channel. gs is nonzero only if S is even. (72? g^i and g^ are known while go has not 
yet [3 El [9l Ho]- However, many features of the condensate depend strongly on g^ 
(say, the phase- diagrams [101 El lEl Il3] and the spin-evolutions [E]). Therefore, the 
determination of g^ is important for a thorough and clear description of this condensate. 




(1) 
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As an application of our idea of repeating collisions in a trap, the ^^Cr atoms have been 
chosen as an example to see to what extent the interaction can be thereby clarified. 

In the beginning, two narrow and deep potentials |ma;^|r±ap are preset at ±a 
(say, two magnetic traps), where a is lying along the positive Z-axis. Each potential 
contains a Cr atom in the ground state of the parabolic confinement. Both atoms 
are polarized but in reverse directions. The upper (lower) atom has spin-component 
II — 3 (—3). Thus the magnetization Ms of the system is zero, and the two atoms 
are localized. Instantly, the two preset potentials are replaced by a broader potential 
^muj^r^ (say, an optical trap) centering at the origin. Then, the previously localized 
atoms begin to evolve. Since Ms = 0, the S = component must be included and will 
be affected by go. Therefore, by observing the evolution, the knowledge on go might be 
extracted. 



2. Hamiltonian and initial state 



Introduce R = (ri + r2)/2 and r = r2 — ri for the cm. and relative motions. Introduce 
?iuj and yjhjmw as the units of energy and length, respectively. The symmetrized and 
normalized initial state 



I + P12 r/2a\3/4 



\2ttJ 



X3(l)X-3(2), (2) 

where P12 denotes an interchange of 1 and 2, a — ooa/ou, 6 is the angle between r and 
the Z-axis, a — |a|, and x^^ is a spin-state of an atom with component ji. Then the 
evolution is governed by the Hamiltonian 

Hevol — Hr + (3) 

where Hr = -|V|-F-R^ and Hr = -V^-F|r^ and Vu = Vs+Vm- The eigenstates of Hr 
and Hr, denoted as 4>nl{R)Ylm{F() and 0ni('")^mj(^)) a-re just the harmonic oscillator 
states. The eigenstates of Hr -\- V12 can be expanded in terms of basis functions as 

^J2^^'rMr)(lS)j, (4) 

7 

where 7 represents the set n, /, and S, {lS)j denotes the coupling of / and S into the 
total angular momentum J, i is just an index of the J-serics. Due to Vdd, I and S are 
not conserved, but J is. Due to the boson statistics, / -|- S must be even. 

It turns out that Vu is rather weak in our case. Consequently, each eigenstate of 
Hr + V12 is close to an eigenstate of Hr- This fact leads to a great reduction of necessary 
basis functions in the expansion. When a set of basis functions has been chosen, the 
associated matrix elements of Hr-\-Vi2 can be derived as shown in the appendix. Carrying 
out the diagonalization, the coefficients Cf^ and the corresponding eigenenergy can 
be obtained. In terms of ^atl and ipf , the initial state can be rewritten as 



N 
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Making use of Eq. (jl]), and equating and ([2]), it is straight forward to obtain 
the coefficients Bj^ and bf . 

3. Time-dependent density and the repeating collisions 

The time-dependent solution for the evolution is 



N 



(6) 

J,i 

where t = ut. 

From Eq. ([6]), all the information on the evolution can be extracted. The main 
feature of the evolution is the occurrence of repeated collisions, as we shall see, and 
the effect of interaction on each collision can be accumulated and therefore enlarged. 
Thereby the strength go can be evaluated. 

We firstly extract the time-dependent density from "^{t) as 



p{r,e,t) = 27ir^ j c;R^*(t)^(t) 



Ji7, J'i'7' 



0^1 L^iiyO^ Oj^ ^v^-Ms; SMs^l, 



yJ,0 

-Ms; SMs^l,-Ms; SMs 

Ms 



■ (l)n'l'ir)(pnlir)\Yi>_Ms\\Yl-Ms\, (7) 

where the Clebsch-Gordan coefficients have been introduced. It satisfies 

dr sine de p{r,e,t). (8) 



When the parameters are given as w = 2000 x 27r, a = 1.5, a = 2, and Qq = —gQ/2, 
the variation of p in the earliest stage of evolution is shown in Fig. [H Since p does not 
depend on the azimuthal angle and p(r, n — 9,t) = p{r, 9, t), it is sufficient to be plotted 
only on a quarter of z-x plane, where z and x are the components of r/2. Since the cm. 
is always distributed close to the origin, we have roughly r2— r/2 and ri ~ —r/2. Thus 
the distribution of an atom can be understood from Fig. [TJ When the evolution begins, 
the two atoms located at opposite ends of the broad potential collide straightly with 
each other (see[T^ and[T]D). When t ^ 7i/2u (Jit), the two atoms keep close to each other 
with a distance ~ 1.6, and they are both distributed around the center. Afterward, the 
atoms begin to separate ([T]i is very similar to [TJ)). When t = tc/u, the profile (not 
yet shown) is very similar to [1^. Thus the first round of head-on collision has been 
completed, and the second round will begin successively. If we remove V12 from Hevoi, 
the factor {E^,' - Ef) in Eq. ([7]) would become an integral multiple of hui, thereby the 
above process would be exactly periodic with the period 27i/u. In fact, the collision as 
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Figure 1. (Color online) p{r,9,t) plotted in the z-x plane, where z and x are the 
components of r/2. t is given at 4 values in the interval (0 to 0.6tt/uj = 0.00015 sec). 
In the panel, the area with the largest p{r, 9, t) is marked by a x . The parameters are 
given as w = 2000 x 27r, a = 1.5, |a| = 2yjh/mw = 311 nm, and 50 — — ff6/2. 



shown in Fig. [T]is essentially determined by H^. In the early stage, V12 causes only a 
very small perturbation. However, as we shall see, the effect of each collision can be 
accumulated when the time goes on. 

The number of collisions that the two atoms have experienced within t is ~ tcj/vr. 
When t = 0.07 sec and u = 2000 x 2tt, this number is ~ 280. When t is close to 0.07 sec, 
p is shown in Fig. [2l At the first glance, |2^ is similar to [1^. However, the peak in [2^ is 
considerably lower than that of [1^ implying that the atoms are not well localized as in 
the beginning. In fact, p has spread widely in[2^ and contains a smooth peak at x = 1.6 
and z = lA (this smooth peak is too low to be seen in the figure). The two atoms are 
closer to each other in [2b than in [lb. The density varies more vigorously along ^ in [2h 
than in [It. [2]d and [2li do not have the approximate similarity as shown previously in 
[Tb and [Hi. From [2|i and [2f, we know that a round of collision has been completed in 
the interval from 0.07 to 0.07 + tt/w. The evolution in this round is explicitly different 
from that in the first round due to the accumulated effect of V12. 

The accuracy of the above numerical results depends on the number of basis 
functions, which is determined by the scopes of n (from to rimax), I (from to /max), 
and N (from to A^max)- When nmax = A'max = 12 and /max = 14, the associated results 
are found to be nearly the same as those by using rimax = Amax = 10 and /max = 12. 
Thus we believe that the former choice is sufficient. 
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Figure 2. (Color online) The same as Fig. [l]but t is given in the interval (0.07 to 
Om + n/uj). 



4. Time-dependent probability of the spin-component and the evaluation 

Due to the spin-dependent interaction, spin-flips will occur during the evolution. From 
Eq. ([H]) the time-dependent probability of the spin-component of an atom in n is 

P,it)= Yl Sn'Jincos[{Ef: -Ef)T] 

Ji7, J'i'7' 

• ^(2A + 1)^(23' + 1){2S + 1) 

X 

■ W{13J'3; XS')W{13J3; XS) 

where both the Clebsch-Gordan and Wigner coefficients [15] have been introduced. To 
show the convergency of numerical calculation, the dependence of -P^(t) on the number 
of basis functions is given in Tab. [B 

It is clear from the table that, when rimax = -?Vmax = 12 and /max = 14 as we have 
adopted, qualitatively accurate results can be obtained. 

Due to the symmetry inherent in \l/7 and in Hevoi, P-tJ.{t) = Pfiit)- The variation of 
Psi^t) against t is shown in Fig. |3l where qq is given at three values. It is recalled that the 
flrst round of collision ends a.tt = n/u = 0.00025 sec. In this very short time, the change 
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Table 1. Pslt) at three values of t (in sec). The parameters are given as a; = 1000 x 2tt, 
a = 3, |a| ~ 2, and go — ~ge/2. The number of basis functions depends on n,nax and 
^max (-^max — ^^max IS assumcd) hstcd in thc first row. 
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Figure 3. (Color online) The evolution of Pz(t)- The solid, dash, and dash-dot-dot 
curves have go = —g&l'^^ 0, and g&/2, respectively. In |3K, oj — 1000 x 27r, a = 3, 
|a| = 2, while in[3b, u = 2000 x 27r, a = 1.5, |a| = 2. 

of P^it) is negligible. However, after hundreds of repeating collisions, the dependence of 
P^it) on qq can be clearly detected as shown in the figure. Therefore, P^{t) can be used 
to evaluate 5^0- In particular, the dependence will become more explicit if u is larger 
(comparing |3b with[3|i). The evolutions of P2{t) and Poif) are shown in Fig. HIi andHb, 
respectively, where is given at seven values. One can see that there is a small peak in 
appearing in the early stage of evolution if qq is negative. Its height depends on how 
negative go is, and it would disappear if is positive. Therefore, the existence or not 
of this peak can be used to judge the sign of g^. Together with the height and location 
of the second peak in Fig. |3] (or the first peak of Hb), the strength g^ can be known. 
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Figure 4. (Color online) The evolution of P2{t) (a) and Po{t) (b). The curves from 
"1" to "7" have go from —ig^/A to Sge/'^, respectively, with an increase 56/4 in each 
step. Each curve has shifted up by 0.005 more than its adjacent lower neighbor to 
guide the eyes. The other parameters are the same as in Fig. [1] 

Since Vdd can alter /, the spatial structure can thereby be altered. On the other 
hand, Vdd can also alter S; therefore spin-evolutions are also affected. We found that, 
if Vdd is reduced (strengthened), the spin- flips would occur less (more) probable. E.g., 
when Vdd is changed to (3Vdd, the values of the first minimum of the solid curve of Fig. [3b 
would be 0.462, 0.448, and 0.298, respectively, if /3 = 0, 1, and 10. It implies that a 
strong dipole force will cause strong spin- flips. This is a notable point. On the other 
hand, recent experimental progress suggests that condensation of molecules with large 
permanent dipole moments, such as OH [161 E], RbCs [18j, KRb |19], and NH |20], 
may be achieved. These systems would have very strong dipole interaction, 10^ or more 
times stronger than in chromium. Therefore, distinguished phenomena of spins caused 
by the very strong Vdd are expected. 

In order to understand better the phenomenon of spin-flips, the momentum-spin 
correlation is studied in the following. Let the coordinates ri and r2 be introduced via 
the Talmi-Moshinsky (T-M) coefficients as 

= %Tinil2^^ihin)(i)n2l2{r2)[Yl,{h) Yi^{r2)]lm, (10) 

nilin2l2 

where Zi and I2 are coupled to I and m, and 4>ni{ri)Yimi{fi) are the eigenstates of 
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Hr- = —^^ri + I'^i- The analytical expression of T-M coefficients has been given 
in the refs. [211 |22l [23] . Inserting Eq. into (Q, we have 

nili,n2l2,'n2im2,fJ.i' 

■ (f>mh (n)0n2/2 (^2)yiimi (^l)yi2m2 (^2) 

■X^(1)X.(2), (11) 

where F^J^ ^^^^ mim2('^) '^^'^ known by comparing Eq. ( ITT]) with (|6]) and (fTO]) . 

Insert Eq. ( ITT]) into the normality (^(r)|^(r)) = 1, it is straight forward to obtain 
the time-dependent probability of finding a particle in Xu state at r2 and r as follows 

«'(r2)</>n2«2(^2) 

712^2 :"2'2, "12 

■>^/>2 (^2)1^^2-2(^2), (12) 

where 

^n'2l'2,n2l2,m2i'^) ^ ^ '^mi+m2, -A'-i' [^n^i/i ,n^Z^,mim2 (''")] 
nilimifi 

ni/i,n2/2i™i"i2*' ' V / 

This probability fulffils Ylu I ^^"i -^^^(^2, t) = 1- 

For observing the momentum-spin correlation, we define 

U{P) = \/f (-^)' / ^^(^^ Ji{pr)^ni{r), (14) 
where the spherical Bessel function has been introduced. Then from Eq. ( I12p , we obtain 

Qu{P,r)= Yl ^ny2."2«2,m2(^)Q«^(P)W2(p) 

n'2l'2,n2l2,m2 

■yilm2iP)yhm2{p), (15) 

which is the probability of finding a particle in z/ at p and r. It fulfills 

Ei./c^p QAp^^) = 1- 

In order to observe Q^, at a given time Tc, the broad trap has to be suddenly 
cancelled at Tc. Then, the particles begin to go out and the successive evolution is an 
expansion and is governed by the Hamiltonian Hjree + Vu, where Hfree = ~|(Vi-|- V^)- 
During the expansion, collisions happen scarcely. We have already seen that the effect 
of V12 within a few collisions is negligible. Therefore, it is safe to neglect V12 during the 
expansion, thus the wave function after Tc is \['(r > Tc) ~ e-'^^freeiT-Tc)^^ (^^^^ ^ \E'j,f.ee(T). 

It is straight forward to see that the weight of \l//j.ee(T) in |pi,P2) is the same as 
that of ^^(tc). Therefore, after the cancellation, the probability Q,, depends only on Tc 
but not on r. Thus, the time of observation r is not relevant, but the choice of is 
essential. 

An example of Qsi^p, Tc) in the earliest stage of evolution is given in Fig. [Sti. Curve 
"1" describes the momentum distribution of the initial h.o. ground state over the Z- 
component of p. When the evolution begins, the distribution shifts to the left as shown 




Figure 5. (Color online) The distribution of Q3(p, Tc) over the Z-component of p. Tc 
is the time to cancel the trap, and is given at five values as Tc — {i — Ijir/A + rx, is 
marked by the associated curve, and tx = (a) and 2007r (b). The unit of momentum 
in this paper is y/mhuj. The parameters are the same as in Fig. [T] 
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by "2". It implies that the component with u = 3 is moving from the initial end 
toward the center. When r = 7r/2, a sharp peak appears at p^ = —2 as shown by 
"3". Meanwhile, as we know from Fig. [TJ:;, the particle is close to the center, and it 
has a maximal momentum pointing to the opposite end. Afterward, the magnitude of 
momentum begins to decrease. When r = vr, the average momentum is again zero as 
shown by "5" which overlaps "1". Then the process proceeds in reverse direction, and 
repeatedly. When r is larger, the periodicity will be gradually spoiled by the interaction. 
This is shown in Fig. |5b, where "5" overlaps "1" no more. Instead, the peak of "5" 
shifts a little to the left. It implies that the peak has not yet completely arrived at the 
opposite end, its arrival is delayed. 

The momentum-spin correlation can be revealed by QuiPj^c)- However, we are 
more interested in a quantity which is easier to be observed. Therefore, we define 
Qi{Op,Tc) = f p"^ dp QuiPyTc), where 9p is the angle between p and the +Z-axis (the 
azimuthal angle of p is irrelevant). This is the probability of finding a particle in u and 
emitting at a specified direction 9p at Tc. It fulfills 27i J sinOp dOp Qi{Op,Tc) = 1. This 
quantity is much easier to be experimentally measured because only the orientation of 
spin is concerned, which can be measured via a Stern-Gerlach device setting towards 
the outgoing angle 9p. 

Ql is sensitive to (therefore both the time of creation and the time of cancellation 
of the broad potential should be precise). For examples, when Tc = 69.57r (meanwhile the 
average positions of both particles are close to the center), against 9p with given 
at three presumed values are plotted in Fig. [6|l, where the momentum is distributed 
around 6p = 0. Thus the particle is rushing back to the initial end disregarding how 
go is (meanwhile the u = —3 particle is rushing to the opposite end). If Tc = TO.Svr, 
the peak of [6^ would appear at the right (not yet plotted). When = TOtt (meanwhile 
the particle is close to the initial end), Qg depends slightly on 9p as shown in [6]d, the 
dependence is not very sensitive to go. 

On the other hand, due to the spin-flips, the |z/| 7^ 3 components are created and 
the creation is sensitive to the interaction. For the case of the u = 2 component plotted 
in int, the emission with a smaller 9p would be more probable if go is negative, but less 
probable if 5^0 is positive. The probabihty Qi{0, GQ.Svr) would decrease by 7 times if go 
is changed from —ge/4: to 5'6/4. Thus, information on the interaction can be extracted. 

When Tc increases, Q2 would also increase gradually. This is shown in [Gt, [6^, and 
Ef, where Tc = 2k'K — 7r/2 (note that the ordinates of these figures have different scales). 
An example with = 2kTi is shown in[6ll. Comparing [6li withlHt, one can see how Q2 
is changed after the interval ti/2. 

The features of other |z/| 7^ 3 components are more or less similar to those oi v = 2. 
Information on the interaction can also be extracted from them. All the curves in Fig. [6] 
will undergo a left-right reflection if Tc is changed from 2fc7r — 7r/2 to 2kTT + 7r/2, or if u 
is changed to —v. The reflection is nearly exact. 

It is emphasized that the dependence on interaction is hardly found in the early 
stage. For an example, when Tc = 27r (meanwhile both particles nearly completely 
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Figure 6. (Color online) Angular distribution of Qi{Op,Tc) against 9p. ~ 3 for (a) 
and (b), and v = 2 for the others, is fixed and marked on the panels, is given 
at three values —56/4, and ge/i, respectively, and the associated curves are in solid, 
dash, and dash-dot lines. The other parameters are the same as in Fig. [T] 
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return to their initial positions after two rounds of collisions) , the curves of will 
become a set of horizontal lines and very close to each other, and therefore no valuable 
information can be extracted. Thus the introduction of the trap, that leads to repeated 
collisions, is necessary so that the weak effect of interaction can be accumulated. 
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Figure 7. (Color online) Qi/Qi with 9p = tt/2 and Tc = 69.57r against go/ge- v is 
marked by the curves. The other parameters are the same as in Fig. [1] 



When 9p and Tc are both fixed, an example is shown in Fig. [7] where the ratio 
Qi/Qi with 9p = 7r/2 and = GQ.Svr is plotted against go/ge- The curves with z/ and 
—u overlap with each other entirely (however, if 6p ^ 7r/2, they do not). It is clear from 
the figure that g^ can be uniquely determined if a few ratios can be measured (say, those 
with V = 2 and 0). 



5. Conclusions 



In summary, an idea is proposed to study the trapped 2-body scattering. The two 
particles are first localized, then they collide with each other repeatedly in a trap, 
then they escape. Related theoretical derivation and numerical calculation have been 
performed to study the coordinate-spin and momentum-spin correlations. The ^^Cr 
atoms have been chosen as an example. Due to the repeating collisions in the trap, 
the effect of the weak interaction can be accumulated and enlarged, and the strength 
go can be thereby determined. It is expected that the approach might open a new way 
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for studying various spin-dependent interactions among atoms (molecules) with nonzero 
spin, in particular for the studies of very weak forces. 
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Appendix A. Matrix elements of V12 between the basis functions 



{(l)n'i'ir)il'S')j>\Vu\Mr){lS)j) 



5.V 



J', J 



--gs5i'o5io5s's(pn'oiO)(PnoiO) - 252x/5Cdv/(2^' + 1)(2^ + 1)(2/ + 1) 



1 1 2 



Ciii,oClo%,oW{l2jS'-l'S){ 3 3 5V I -<p.Mr)Mr) , (A.l) 

3 3 S' 



dr 



where the Clebsch-Gordan, Wigner and 9-j symbols [15J are introduced. Furthermore, 
due to the Boson statistics, both s + I and s' + I' should be even. 
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